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THE CLUSTER-MERGER SHOCK IN IE 0657-56: FASTER THAN THE SPEEDING BULLET? 
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ABSTRACT 

Shock waves driven in the intergalactic medium during the merging of galaxy clusters have been observed 
in X-ray imaging and spectroscopy. Fluid motions inferred from the shock strength and morphology can be 
compared to the cold dark matter (CDM) distribution inferred from gravitational lensing. A detailed reconstruc- 
tion of the CDM kinematics, however, must take into account the nontrivial response of the fluid intracluster 
medium to the collisionless CDM motions. We have carried out two-dimensional simulations of gas dynamics 
in cluster collisions. We analyze the relative motion of the clusters, the bow shock wave, and the contact dis- 
continuity and relate these to X-ray data. We focus on the "bullet cluster," IE 0657-56, a near head-on collision 
of unequal-mass clusters, for which the gas density and temperature jumps across the prominent bow shock 
imply a high shock velocity 4,700 km s -1 . The velocity of the fluid shock has been widely interpreted as the 
relative velocity of the CDM components. This need not be the case, however. An illustrative simulation finds 
that the present relative velocity of the CDM halos is ~ 16% lower than that of the shock. While this conclusion 
is sensitive to the detailed initial mass and gas density profile of the colliding clusters, such a decrease of the 
inferred halo relative velocity would increase the likelihood of finding IE 0657-56 in a ACDM universe. 
Subject headings: dark matter — galaxies: clusters: general — galaxies: clusters: individual (IE 0657-56) — 
intergalactic medium — shock waves — X-rays: galaxies: clusters 



1. INTRODUCTION 

During the growth of structure in the universe, nonlinear 
cold dark matter (CDM) halos collide and merge hierarchi- 
cally. When halos collide, the dark matter of one can, if 
collisionless, pass through that of the other, unlike the fluid 
baryonic gas in the halos. Since the collisions are frequently 
supersonic for this gas, it passes through shocks where it 
heats and virializes in the final halo. Currently, an X-ray 
surface brightness and temperature discontinuity identifiable 
as a merger shock with Mach number significantly above 
unity has bee n detected with Chandra in two galaxy cl usters: 
IE 06 57-56 jMarkevitch et al]l200l lMarkevitcr]l2006l) and 
A520 (Ma rkevitch et alJ l2005). In both cases, X-ray maps ex- 
hibit a bow-shock like temperature and density jump. For IE 
0657-56, the spatial segregation of the X-ray emitting plasma 
from the peaks of the mass distribution detected with gravita- 
tional lensing has been interp reted as the first d irect proof of 
the existence of dark matter dClowe et alj|2006l) . 

The strength and geometry of these shock fronts, and the 
structure of the contact discontinuity, depend on the kinetic 
energy and the detailed gravitational and gas density profiles 
of the merging components. If one could recover the colli- 
sionless CDM dynamics from the fluid dynamics of the X-ray 
emitting gas, then the kinematics of merger shocks would be 
a powerful probe of the nature and cluste ring of CDM (e.g., 
iMarkevitch et alJl2005l:IClowe et alJl2006h . However, the two 
components are coupled only by gravity and behave differ- 
ently in the merging process. Recently, it has been assumed 
that the gas shock velocity inferred from the X-ray measure- 
ments of IE 0657-56 equals that of the colliding clusters, im- 
plying a rare merger event that might be in conflict with the 
statistical expectations of the CDM model. We critically ex- 
amine this assumption by simulating gas dynamics in clus- 
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ter mergers and studying the response of the intergalactic 
medium (IGM). 

In the cluster-merger in IE 0657-56 1 , transverse motion of 
the merger shock (w 4,700 km s -1 ; Markevitch 2006) greatly 
exceeds relative radial motion of galaxies in the two compo- 
nents (w 600 km s" 1 : iBarrena et alj|2002l) . This suggests that 
the shock was driven in the IGM of the larger component by 
a near head-on supersonic passage of the smaller component. 
Inside the bow shock, another bow-shaped discontinuity with 
a reverse temperature jump (lower temperature on the con- 
vex side) is seen; this has been interpreted as the contact dis- 
continuity separating the shocked an d ram-pressure stripped 
IGM of the two merging components. Takizawa (2005, [20061) 
recovered the three basic features (a bow shock, a contact 
discontinuity, and ram-pressure stripping) in hydrodynamical 
simulations of cluster collisions, but did not address the rela- 
tive kinematics of the CDM and the fluid component. 

We study the formation and propagation of the merger 
shock wave and the evolution of the contact discontinuity, and 
relate these observable features of X-ray surface brightness 
and temperature maps to CDM kinematics. In § [2] we de- 
scribe the numerical method, report the general features of the 
simulations, and analyze a particular run that approximates 
the observed IE 0657-56 (the conclusions also apply to other 
cluster mergers). In § [3] we discuss implications for CDM 
clustering. T he standard co smological parameters consistent 
with WMAP ( Sperg el et ai1l20( )6) are assumed throughout. 

2. SIMULATIONS OF CLUSTER MERGER SHOCKS 

2. 1 . The Algorithm and Initial Conditions 

The simulations were carri ed out with the tho roughly-tested 
Eulerian code ASC FLASH dFryxell et al.ll2000h in two spa- 
tial dimensions. Its adaptive mesh refinement capability al- 



1 With an average X-ra y temperature of 14. 1 ± . 2 ke V I Markevitclj^etalJ 
1 200a IMarkevitcr] 12006: Andersson et al. 2006; Markevit cT& Vikhlininl 
120071) . the cluster IE 0657-56 at redshift z = 0.296 fTucker et al]|1998F "is 
the hottest and X-ray brightest known cluster. 
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lowed us to simulate a large spatial domain with dimensions 
of 20 Mpc and simultaneously achieve a spatial resolution of 
~ 1 kpc at fluid discontinuities. The simulations were carried 
out in cylindrical (r,z) coordinates in an inertial frame. Axial 
symmetry restricts the simulations to head-on cluster colli- 
sions. To avoid artifacts which can result when a perfectly 
head-on collision occurs with cuspy halo profiles and to ap- 
proximate the conditions in non-head-on collisions, the cen- 
tral dark matter and gas density of the larger cluster was artifi- 
cially reduced. The gravitational potential of each of the clus- 
ters is spherically symmetric and a fixed function of distance 
from the center. The clusters are placed at separation D12 and 
are allowed to move under each other's gravity. The accelera- 
tion of each cluster is evaluated at its center. The effect of the 
Hubble flow on cluster motion is ignored because it only af- 
fects the early infall, while we are interested in the dynamics 
of the two halos and their gas content at pericenter passage. 
We also ignore the dynamical friction drag that reduces the 
halo velocity during pericenter passage. The fluid flow and its 
relation to CDM dynamics within a few hundred kiloparsecs 
from the center of the smaller halo should not be different in 
the presence of dynamical friction, because the gravitational 
field of the dynamical wake that trails the smaller CDM halo 
will not have strong gradients (tides) on such scales. Finally, 
we ignore the tidal stripping of the smaller CDM halo which 
could be substantial on larger scales than are relevant here. 

For the mass distribution a ssociated with the gravita- 
tional potential we consider the iNavarro et al.l (1 19971) profile 
("NFW") pnfw = po{r/r % T l {\ +r/r s y 2 , and the "core" profile 
Pcore = />o(l +r/r t ) ~ 3 . The profiles were normalized such that 
the mass enclosed within a fixed radius r^w equals the mass 
of the cluster M500. The gas density profile was allowed to 
differ from that of the CDM profile. At the start of each sim- 
ulation, monatomic gas with adiabatic polytropic equation of 
state (7 = 5/3) is set in hydrostatic equilibrium within each 
halo. Explicit viscosity, heat conduction, cooling, and self- 
gravity of the gas are ignored. Runs were carried out for a 
variety of merger velocities and density profiles that included 
NFW and core CDM profiles, as well as cuspy (divergent) 
and non-cuspy gas density profiles. In divergent CDM density 
profiles, a small constant-density core (20 kpc) in hydrostatic 
equilibrium was applied in the center to resolve the fluid gra- 
dients with multiple resolution elements. The cluster gas frac- 
tion was /gas^b/^m, where fib/O m « 0.16 is the cosmologi- 
cal baryon abundance and / gas ss 0.7 ( Afsh ordi et al.l l2006). 2 
The mean molecular weight was /1 = 0.59, corresponding to 
75% hydrogen and 25% helium by mass, fully ionized. Nu- 
merical convergence was ascertained by resimulation at in- 
creased spatial resolution. 

2.2. General Features of Cluster Merger Shocks 

Two prominent features seen in all simulations are a bow 
shock wave, and a contact discontinuity (cold front). The for- 
ma tion of a shock w a ve d uring cluster merger was discussed 
by iMarkevitch et al.l (l2000h ; our results are consistent with 
their picture. Opening angle and the radius of curvature of the 
nose of the bow shock are sensitive to the details of the simu- 
lation, but both are larger than those expected in steady-state 
bow shocks driven by spherical hard sphere s moving with a 
constant velocity in a uniform medium (e.g.. lFarris & Russelll 

2 All gas-dynamical results are independent of this assumed value of / gas , 
except that the X-ray surface brightness oc /g as . 
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FIG. 1. — Temperature in the meridional plane after cluster pericenter 
passage. The larger CDM halo is located at (r, z) = (0,0). The circles indicate 
the scale radii r s of the two halos. 

[19911 IVikhlininetaIll200ll and refer ences therein). Mate- 
rial on the convex side of the contact discontinuity originates 
in the smaller cluster. During the collision, a reverse shock 
propagates into the smaller cluster; a fraction of the smaller 
cluster traversed by the reverse shock is heated, but the dense 
core of the cluster remains cold as the reverse shock weakens 
inside it. The wings of the contact discontinuity are Kelvin- 
Helmholtz (KH) unstable (in real clusters, the KH instability 
will be modified by magnetic stresses). A nonlinear instabil- 
ity disrupts the coherence of the surface of the discontinuity 
and creates a narrowing in the axial diameter of the disconti- 
nuity (a neck) just behind the nose, giving it a mushroom-like 
or bullet-like appearance, as evident in the temperature map 
in Figure Q] At later times (not shown), the Rayleigh-Taylor 
instability e ats into and disrupts the bullet from the rear, as 
predicted by IMarkevitch & Vikhlininl d2007l) . 

In most simulations, during pericenter passage, the shock 
and contact discontinuity are slower than the relative veloc- 
ity of the two CDM halos (shock velocity is calculated in the 
local rest frame on the upstream side of the shock). This is 
due to the ram-pressure force acting on the gas in the center 
of the smaller cluster but not on the halo. In many runs some 
degree of ram-pressure stripping is evident in the center of 
the smaller cluster (we call the smaller cluster ram-pressure 
stripped when cold gas originating in the smaller cluster has 
been expelled from the center of the smaller halo). However, 
later, as the halos climb out of each other's gravitational po- 
tential well and decelerate, but the shock and the contact dis- 
continuity do not decelerate appreciably over a longer period 
(Fig.lJJ;). The lack of deceleration of the shock and the cold 
front was noticed, though not quantified, in previous s imula- 
tions of cluster mergers (Ma rkevitch & Vikh linin 2003, and 
references therein). It could result from the drop in ram pres- 
sure as the cold bullet propagates into a thinning cluster at- 
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FIG. 2. — (a) Velocity of the shock (solid) and of the contact discontinuity 
(dashed) relative to the pre-shock gas upstream, and relative velocity of the 
CDM halos (dot-dashed), all as functions of the time-varying separation be- 
tween the halos D\i = zi — Zi . (b) The density (solid) and temperature jump 
(dashed) at the shock, and the shock Mach number (dot-dashed), (c) Distance 
of the shock (solid) and the contact discontinuity (dashed) from the center of 
the smaller halo. 

mosphere and from the gravitational tide of the halo. After 
the cold gas originally in the smaller cluster has been ram- 
pressure stripped, the stripped gas lags behind the smaller 
halo, but eventually its velocity becomes larger than that of 
the halo, which implies that the gas can later catch up with 
the halo again. 

Before proceeding to discuss a specific simulation approxi- 
mating the observed bullet cluster, we strongly emphasize that 
the shock kinematics and morphology are extraordinarily sen- 
sitive to the parameters of the simulation. In particular, a small 
variation (~ 25%) in the cluster mass ratio or density profile 
can lead to vastly different degrees of ram-pressure stripping, 
or the opening angle of the bow shock. 

2.3. A Model of IE 0657-56 

General results discussed in § 12.21 are illustrated by a spe- 
cific run that came the closest to reproducing the observed 
properties of IE 0657-56. Cluster masses in the run were 
M 500 ,i = 1.27 x 10 15 M o and M 500 ,2 = 2.54 x 1O 14 M , cluster 
radii were ^500,1 = 1.5 Mpc andM5oo.2 = 1 Mpc, and the cluster 
scale radii were r Sj i = 500 kpc and r s 2 = 333 kpc. 3 The density 

3 The mass ratio of the two halos in the simulation, M sgoi/Msgog = 5, is 
about a half of the ratio ~ 10 inferred from lensing data I Clowe et aT] |2004l 
2006; Bradac et al. 2006). Some fraction of the smaller cluster's initial mass 
will have been tidally stripped during the merger, possibly reconciling our 
initial conditions with lensing. The total luminosity L lot of galaxies in each 
cluster can be u s ed to estimate the cluster mass via the relation Lt ot <x M^ Q0 
)Lin et all 120041 IMiller et all 120051 ; ICoorav & Milosavlievic 2005), where 
j3 < 1 is a bandpass-dependent factor. For the S-band we adopt (3r = 0.65. 
Barrena et al. 12002) report L lot .i r = 1O 12 L and L,ot,2,ff = 0.2 X 10 12 L o , im- 
plying a mass ratio of A/500,1 /^500 2 ~ 12, in agreement with lensing. The 
difference of the mass ratios may also be attributed to a slightly non-head- 
on collision in IE 0657-56. Indeed, the lensing maps are consistent with a 
northerly passage of the smaller cluster inducing the observed SW— NE tidal 
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FIG. 3. — Temperature profile on the central axis (left panel, solid line); 
emissivity- weighted integrated temperature profile (left panel, dotted line); 
and 0.8-4 keV surface brightness on the central axis (right panel). The cen- 
ter of the larger CDM halo is located at z = 0. The bow shock at z ~ 0.94 Mpc 
is split into a main shock and a weaker subshock. 

profile of the larger cluster was the non-divergent p core with 
central density « 2.3 x 10~ 3 M Q pc~ 3 (again, to approximate 
the conditions in a non-head-on merger with pericenter pas- 
sage at a distance ~ r s .1), and, in the smaller cluster, it was the 
NFW profile pnfw- The gas density in both clusters was pro- 
portional to the CDM density, /9 gas = (1 + ^m//gas^b)~ I /OCDM- 
The halos initiated infall from relative rest at separation D12 = 
4.6 Mpc and their maximum relative velocity at pericenter 
was y max » 5,270 km s -1 . 4 

Figure Q] shows a snapshot of the meridional temperature 
distribution when the distance of the shock to the larger CDM 
halo is about 1 Mpc. The larger and smaller cluster are mov- 
ing vertically upward and downward in the figure, respec- 
tively. The cold gas bullet is ~ 100 kpc in radius, excluding 
visible fingers of cold gas embedded in the hot post-shock 
medium. The center of the smaller halo is located ahead 
(below) the contact discontinuity, and almost exactly at the 
location of the bow s hock, as found in X-ray and lensing 
maps of IE 0657-56 dClowe et al.ll2004l 120061 iBradac et all 
2006). This is also evident in Figure |2j; showing the position 
of the bow shock and the contact discontinuity relative to the 
position of the smaller halo. Ram-pressure stripping at the 
center of the smaller cluster occurs at D12 = 0, but the bow 
shock and the contact discontinuity catch up with the halo at 
D[2 « 1.2 Mpc and w 1.5 Mpc respectively. 

The pre- and post-shock axial temperatures are 12 keV and 
34 keV, respectively (see Figure[3]l. Figure|2j? shows the evo- 
lution of the density and temperature jumps across the shock, 
and its Mach number, as a function of the halo separation. 
The Mach number is rs 3 when the separation equals the pro- 
j ected separation of lensing density peaks D12 ~ 720 ± 25 kpc 
(ICloweet all [2004. 2006; Bradac etal J I2006T) . The temper- 
atures and shock s trength are compatible with the results of 
iMarke vitch (2006), who, for a 7 = 5/3 polytropic equation of 
state and instant electron-ion equilibration, derived the Mach 
number from deprojected density and temperature jumps in 
the bow shock in IE 0657-56. The separation of the shock 
and the contact discontinuity is ~ 70 kpc, somewhat smaller 
than in the observed system. Figure [4] shows the predicted 
X-ray surface brightness map for the simulated cluster. 

Figure^ shows the evolution of the velocities of the shock, 
the contact discontinuity, and the relative velocity of the two 
CDM halos. While the cluster velocity decreases as the small 

distortion in the larger cluster. Then, the lower density column met by the 
small cluster would result in equivalent shock dynamics if the smaller cluster 
were smaller than in the simulation. 

4 All dimensionless numbers (such as the shock Mach number, and the gas- 
to-dark matter velocity ratio) and distances (such as the radius of curvature 
of the shock) are invariant under a uniform scaling of the CDM density that 
preserves lengths. The scaling preserves the mass to temperature ratio. 
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FIG. 4. — Surface brightness map of the simulated bullet cluster in the 0.8 - 
4 keV band, assuming a frequency-independent effective area of 400 cm . 

halo climbs the potential of the larger halo, the shock veloc- 
ity remains constant. The lack of deceleration of the shock 
can be explained by a combination of factors; in addition to 
gravitational forces, the negative density and temperature gra- 
dients into which the shock propagates tend to increase its 
strength. At D12 = 720 kpc, the velocity of the shock is V S h « 
4,800 km s" 1 , consistent with V sh ~ 4,740+^" km s~' inferred 
for IE 0657-56 by Markevitch (as cited in iFarrar & Rosenl 
2006, private communication). However, the relative ve- 
locity of the halos in our simulation is much less, V su b ~ 
4,050 km s" 1 . 5 

3. DISCUSSION AND CONCLUSIONS 



Rec ent analyses of I E 0657-56 (e.g., lH avashi & White 
2006t IClowe et al.ll2006t iBradac et all 120061: IFarrar & Rosenl 
2006) assume that the shock velocity equals the instan- 



taneous relative velocity of the CDM halos, and notice 
that it is much larger than t he vir ial velocity of the larger 
cluster. lHavashi & White! d2006l hereafter HW06) and 
IFarrar & Rosenl (120061 hereafter FR06) estimate the proba- 
bility that a halo in the ACDM universe has a subhalo with 
velocity larger than that of the bow shock in IE 0657-56. 
We assess a correction that should be applied to the inferred 
probabilities to account for the relative motion of the shock 
and the subhalo. HW06 calculate the likelihood of find- 
ing a halo-subhalo pair with such relative velocity by ex- 
trapolating results of the Millenium cosmological /V-body 



simulation. They find that it decreases rapidly with V sl ,b, 
log[M(> V su b)/AW s ] = -(Vsub/vioV20o) Q , where vio « 1.55 
and a ks 3.3. The probability depends on the assumed val- 
ues of V sub and 7200. in which HW06 and FR06 differ. 6 
HW06 adopt a shock velocity of V sh w 4,500!^° km s" 1 
( Mark evitch et al.1 120041) as the relative velocity V su b of the 
two halos and adopt V200 ~ 2, 380 km s -1 for the virial ve- 
locity. With these, they infer that the probability of most 
massive subhalo having such a velocity is N\(> V S ub) ~ 1%, 
whereas when they adopt a lower shock velocity consistent 
with uncertainties of V su b = 3,700 km s -1 , the probability is 
Ni(> Vsub) ~ 10%. In our simulation (§ |2~3]l, we find that, 
when the halos are separated by the observed projected dis- 
tance D12 = 720 kpc, the velocity of the shock and that of the 
subhalo differ by ( V s h - V S ub)/ Vsub ~ 16%. Thus, the prob- 
abilities reported by HW06 should be corrected and should 
read « 8% and rs 27%, respectively. FR06 adopt a some- 
what higher shock velocity V s h ~ 4,740 km s -1 and signifi- 
cantly lower virial velocity V200 ~ 1 , 740 km s -1 . These values 
imply a very low probability of w 3.7 x 10~ 7 . 7 The probabil- 
ity is significantly higher but remains very low, « 2.4 x 10 -4 , 
when the 16% correction for the difference between the shock 
velocity and CDM halo velocity is taken into account. 

In conclusion, our 2D hydrodynamical simulations of merg- 
ing galaxy clusters show that, in a case like the "bullet cluster" 
IE 0657-56, the halo collision velocity need not be the same 
as the intergalactic gas shock velocity. While the kinematics 
of the shock is sensitive to the details of the cluster structure, 
the instantaneous shock velocity can exceed the relative ve- 
locity of CDM halos by at least ~ 16%. Any attempt to relate 
shock kinematics to the details of hierarchical clustering and 
halo assembly must take into account the nontrivial response 
of the IGM. Published estimates of the likelihood of finding 
a configuration resembling IE 0657-56 in a ACDM universe 
may require upward revision by a factor of 2 to 700. 
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If the true halo separation is larger than the projected separation 720 kpc, 
the predicted subcluster velocity V su (, will be < 4,050 km s -1 . 

6 We do not test the accuracy of the cited values of V200. but the lensing 
data seem to indicate that HW06's value of V200 is an overestimate. 



7 HW06 and FR06 also report X 5 lower probabilities after counting only 
post-collision systems. 
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